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ABSTRACT 


The Naval Oceanographic Office developed the Integrated Carrier 
Acoustic Prediction System (ICAPS) in order to perform on-scene 
acoustic forecasts for the Fleet. To evaluate the system, bathy- 
thermographic information was obtained for points in the North 
Atlantic during the naval exercise "SEACONEX." Fleet Numerical 
Weather Central (FNWC) provided climatological and synoptic tempera- 
ture profiles along with related acoustical forecasts for the points 
under study. Similarly, ICAPS supplied a combined climatological 
and synoptic temperature profile and related acoustic forecasts for 
each of the points. Comparisons of the two systems were made 
against an actual mean temperature profile and the acoustic fore- 
cast that it generated. It appears that the climatological acoustic 
forecasts were not reliable tactically. Yet, the synoptic forecasts 
by FNWC did yield reliable low frequency acoustic forecasts within 
the limits of computer roundoff error. ICAPS capability for only 
aoe eemmera eure sprotile input proved to yield reliable acoustic 
forecasts, and should be uSed as the real-time acoustic forecast 


system for the Fleet. 
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I. INTRODUCTION 


/,  SeeORBIECTIVES 

The primary objective of this study is to compare the ocean 
thermal structure analyses and acoustical products of Fleet 
Numerical Weather Central (FNWC) with the Integrated Acoustic Pre- 
diction System (ICAPS) thermal structure inputs and acoustic fore- 
casts. This comparison is intended to evaluate the overall 
credibility of both systems. The data used in this study were 
collected aboard the USS GUAM (LPH-9) during the naval exercise 
"SEACONEX" held in the North Atlantic Ocean in June 1973. 

i Orden tO dceconplish) this obyeccive (basic plan including 
discussion of FNWC models, etc.) this paper presents a brief dis- 
cussion of the merits of the climatology used by FNWC and the new 
interim climatology developed by the Naval Sn AO Office 
for the ICAPS system. Additionally the FNWC synoptic temperature 
profile and its related acoustic field are compared to those 
generated by using ee mean bathythermograph (BT) profiles 
obtained at sea. The value of a Single BT in the real ocean is 
discussed. e the case of a single BT profile obtained aboard 
ship, a theory concerning the errors associated with using such 
a profile is advanced. 

A second objective of this study is to appraise the operational 
techniques required for ICAPS. By a detailed discussion of the 


Strengths and weaknesses of ICAPS temporary profile inputs combined 





with suggestions on how to capitalize on the strengths and compen- 
sate for the weaknesses, a pattern of operation will be projected. 
Planning and operating procedures for the use of ICAPS are proposed. 
MPissione are aimed at the procurement and processing of oceano- 
graphic data rather than the tactical utilization of such informa- 


tiOn~es The goal of this portion of the study is to improve the 


effectiveness of the shipboard acoustical forecaster. 


B. ACOUSTIC FORECASTING MODELS 

The Fleet Numerical Weather Central (FNWC) in Monterey, Cali- 
fornia, currently produces and distributes to the fleet two primary 
acoustic forecasting products, the Ship, Helicopter Acoustic Range 
Prediction System (SHARPS II) and the Acoustic Sensor Range Predic- 
tion (ASRAP III) (Navweaservcom, 1971la,b). 

1. Ship, Helicopter Acoustic Range Prediction System (SHARPS IT) 

The Ship, Helicopter Acoustic Prediction System (SHARPS) is 
a computerized sonar range prediction scheme. It 1s designed to 
take environmental analyses and convert them into acoustic range 
predictions. Inputs. to the program include average equipment para- 
meters, platform speed, sonar operating mode, Sensor and target 
depth. Range predictions are based upon a Single-ping, 50 percent 
probability of target detection by an unalerted sonar operator. 

The key input to the SHARPS II program is the vertical tem- 
perature profile. This, along with the climatological salinity 
profile, is transformed by the use of Wilson's (1960) equation 
Peeweoundspecd IntO a Sound velocity profile (SVP). In turn, the 


SVP is used to generate a raytrace pattern considering direct path, 





purely refracted, and the surface and bottom reflected propagation 
paths. The change in sound intensity with range is computed con- 
sidering losses due to spreading, scattering, reverberation, absorp- 
tion, and surface and bottom reflections. Average performance 
characteristics are used for each sonar type, and an estimate of 
the sound intensity as a function of range is achieved. Considering 
that for active sonars a two-way travel distance (distance from 
source to target and return)is needed; the amount of propagation 
loss is double that of the one-way travel. 

SHARPS II claims increased accuracy over previous models. 
This accuracy is attributed to (1) use of the modified Atlantic 
Meteorological and Oceanographic Study (AMOS) equations for surface 
channel, direct path, and sub-surface duct computations; (2) inclu- 
sion of the physics of the operational raytrace and propagation 
loss model techniques for bottom bounce and convergence zone compu- 
tations; (3) inclusion of the basic reverberation technique used in 
the Navy Interim Surface Ship Sonar Prediction Model (NISSM) 
developed by Mobile Sonar Technology (MOST) Committee; and (4) in- 
creased input data = the model including a description of the SVP 
from the surface to the bottom, a detailed sonar system description 
Or <a 1 Pe niente, and area scattering strengths for reverberation 
computations. 

Each prediction is made for a predetermined acoustical area 
in which several assumptions are made: (1) water depth is constant; 
(2) SVP is horizontally constant; and (3) bottom composition and 


Peug@inesssare CONStant. Predictions are made on the basis of 12 
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hourly environmental temperature analyses by FNWC. The reliability 
and variability of that temperature input directly affects the 
reliability and variability ofthe acoustic forecast. 

2. Acoustic Sensor Range Prediction (ASRAP IIT) 

The ASRAP III program is directed toward the requirements 
of the airborne acoustic Sensors, primarily the passive sonobuoys. 
Normally two distinct forecasts are available: a forecast of passive 
propagation loss at discrete frequencies, and a forecast for the 
active airborne sensors, primarily the active sonobuoys. 

ASRAP III uses a sophisticated propagation loss model which 
employs an analytical solution to the wave equation to determine 
the acoustic intensity for various ray paths. The ocean is divided 
into acoustically homogeneous areas in which factors that affect 
sound propagation remain constant in the horizontal. [In general, 
the intensity level at a receiver is a combination of sound arriv- 
ing from the same four paths mentioned in the SHARPS discussion, 
Losses of energy by reflection from a boundary, spreading, and 
scattering are calculated for propagation within the surface duct 
and below it. Of SBtes6, with a passive system only one-way trans-~ 
aon from the source to receiver must be considered. 

The SVP used by the ASRAP program is derived in the same 
manner as that used by SHARPS. Similarly, the reliability of the 
ASRAP acoustical forecast is a function of the reliability of the 


temperature profile. 


aut 





Cae JeGRATED CARREBER ACOUSTIC PREDICTION SYSTEM (ICAPS) 

The Naval Oceanographic Office (NAVOCEANO) and the Naval Under- 
sea Center, San Diego (NUC), have developed an acoustic prediction 
system capable of providing on-scene forecasts. A mobile unit 
capable of rapid forecasts, yet small enough to be used aboard de- 
stroyers, submarines and VP aircraft, this system has added a new 
capability to acoustic prediction. The system, now called the In- 
tegrated Carrier Acoustic Prediction System (ICAPS), is composed 
of a Data General Corporation Nova 800 computer, a Pertec tape 
recorder, a Tektronix 4002A graphic computer terminal with CRT, 
Guaeetektronix hard copy unit (Consultec, 1972). 

In order to ensure that the ICAPS forecast is acoustically valid, 
the system has been programmed with the SHARPS II and ASRAP III 
models. This gives the unit the advantages of the sophisticated 
acoustic models used by FNWC. A significant difference does exist 
between this on-scene approach and that used by FNWC. FNWC uses 
avallable BT reports from the preceding 72 hours and depends 
heavily upon climatology, even in the near surface region. ICAPS 
uses recent Br aera. or composites of them, in the near surface 
region. NAVOCEANO has develcped a preliminary climatology for 
ICAPS which provides the thermal structure at depths greater than 


range of the BT. 


D. TEMPERATURE PROFILE DERIVATION 
1. Fleet Numerical Weather Central 
FNWC computes the thermal structure in the upper 1200 feet 


in the Northern Hemisphere every 12 hours. This analysis is 
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composed of three parts: a sea surface temperature analysis, a 
mixed layer depth analysis, and a sub-surface temperature analysis. 
a. Sea Surface Temperature Analysis 
The present scheme used by FNWC for sea surface tempera- 
ture (SST) analysis is the Fields by Information Blending (FIB) 
technique (Holl,et al, 1971). FIB performs numerical analysis of 
temperature fields by the assimilation and the blending of informa- 
tion from observations and climatological averages of sea surface 
temperatures. It combines information from three stages in time-- 
current, near past and historical temperatures. Details of the 
analysis procedure are given in Appendix A. 
b. Mixed Layer Depth Analysis 
The mixed layer depth (MLD) is computed by the POTMLD 
program, which takes into account hoth climatological and synoptic 
information (Hesse, 1973). Details are given in Appendix B. 
c. Sub-surface Temperature Analysis 
The sub-surface analysis model as described in the new 
FNWC Computer Products Manual (Hesse, 1973) is constructed on the 
premise that the mee of BT's available in any one day is small, 
much too small to make a truly synoptic analysis possible. Thus, 
it seeks to ae use of the BT information that is available, but 
with a strong dependence upon climatology and the most recent SST 
analysis. Details are given in Appendix C. 
ad. FNWC Climatology 
Stevens (FNWC, personal communication) indicated that 


FNWC developed its Atlantic climatology based upon National Oceanic 


Ee 





Data Center (NODC) station data, various atlases and from 500,000 

- 1,000,000 mechanical BT obServations. This data waS weighted so 
that the effects of a small amount of data from any time period 
would not be treated with the same weight as a larger amount of 
data for another time period or area. This step waS necessary to 
ensure that the resulting climatology profile would be representa- 
tive of the entire time period to which it pertained, and would not 


be dominated by any limited portion of the data. 


2. $ICAPS TenperaGure Profile 


Fox (NAVOCEANO, personal communication) stated that either 
a Single BT observation or a hand~-prepared composite may be used 
aS an input. This is merged at 1200 feet by the computer with the 
ICAPS climatology. There are preSently no capabilities for storage 
Gfebr reports for the previous 72 hours nor is there a function for 
the preparation of composite BT profiles in time or space. Also, 
no gross error check (GEC) is performed by the computer. 

a. ICAPS Climatology 

The Atlantic climatology used by ICAPS was developed 

by NAVOCEANO as an interim climatology until a more sophisticated 
model could be developed. Temperature data were taken from NODC 
files and the NAVOCEANO oceanographic atlas. . All available infor- 
mation within each five sogtbe, hak of ocean area was averaged 


at Standard levels to obtain the climatological profile. An 


average MLD for the area was also included (Fox, personal communication). 
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If. COMPARISON OF TEMPERATURE AND ACOUSTIC PROFILES 


Pee bAolS FOR COMPARISONS 

The purpose of this portion of the paper is to compare the tem- 
perature profiles used by FNWC with those used by the ICAPS program 
for the generation of their respective acoustical forecast canes 
discussed in Chapter I, the temperature input to ICAPS normally con- 
Sists of the latest available BT to a depth of 1200 feet. Below 
1200 feet the BT trace is merged by the computer with the ICAPS 
climatology (Consultec, 1972). For the upper 1200 feet FNWC uses the 
available SST and BT reports, blending them in the horizontal with 
climatology and the analysis for the previous 12 hours. The analyzed 
temperature profile is then merged with the FNWC climatology below 
1200 feet (Hesse, 1973). When BT reports are not available, both 
FNWC and the ICAPS program can generate acoustical forecasts using 
only their respective climatologies. 

The data required to compare the FNWC and ICAPS temperature pro- 
Files includes Bee cene nT reports, the FNWC analyzed thermal struc — 
nies, and the FNWC and ICAPS climatologies. The at-sea data were 
Meiected on board USS GUAM (LPH-9) in the North Atlantic Ocean 
during June 1973. These profiles were obtained by the ships and 
helicopters participating in exercise "“SEACONEX." The BT profiles 
were simultaneously used as inputs to the ICAPS program and sent to 
FNWC as inputs to their thermal structure analysis program. The 


data collection points are indicated in Figure 1. It should be noted 


es 
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that due to the nature of the exercise, considerably more data 
feeencolhlected at points By. C, and Dp than at the other points. 
While on board USS GUAM the ICAPS climatology profiles were also 
obtained. Additionally, @Quring June the FNWC thermal structure 
analysis and climatological profiles were collected and held at 
FNWC. All acoustic forecasts were based upon the same salinity 
profile; therefore, that parameter does not enter this study as 
a variable. 

1. Limitations on Data Used in This Study 

Admittedly, because of the severe limitations of the study 

with regards to time (limited to June) and space (given points only), 
conclusions reached must be treated as applicable only as first 
approximations. Another constraint which was present during this 
exercise is that a transiting exercise limits the number of BT's 
that can be taken in an area. Thus any statistical products could 
be in errcer on the same order of magnitude as the expected varia- 
bility in that area. Were all sampies to be taken at a time and 
place experiencing a maximum in deviation and in the same sense, 
then a biased statistical base would be developed. Thus only 
definite trends in this study will be regarded as indicative. Extra- 


polations to other areas or time periods may not be valid. 


B. CLIMATOLOGY COMPARISON 
1. Temperature Profile 
As mentioned in Appendix A, at FNWC temperature data are 
analyzed and merged at designated depths. Since most BT's obtained 


today are valid at least to a depth of 1200 feet, all comparisons 


ay, 





are limited to these levels. The basis for comparison is the simple 
average of all BT's taken at a position during a specific time period; 
this is termed the actual mean temperature profile. Thus, this com- 
parison is restricted to comparing FNWC climatologies and ICAPS ae 
matology for the upper 1200 feet. As differences between the 
climatologies and the actual mean will be studied, the appreciable 
differences will be assumed to occur only in this surface area. 

Since this study deals with deviations from a mean tempera- 
ture profile, it 1S advantageous to separate the geographical areas 
from which data were obtained into areas where tlie vertical tempera- 
ture profiles exhibit relatively small or large standard deviations, 
respectively. ‘This will permit an investigation of any relationship 
between expected standard deviation for an area and the thermal and 
ecoustic reliability of forecasts for that area. 

To properly separate the positions into regions of high or 
low standard deviation, tables of expected standard deviation con- 
structed by FNWC (Lewit, 1972) are used. These tables cover two 
degree squares over all ocean areas. Table I lists the values for 
the positions used in this study. 

Normalization of the standard deviation for each depth is 
accomplished by dividing the standard deviation at that depth by 
the maximum standard deviation of the column. These results are 
plotted in Figure 2 with those areas judged as having large standard 
deviations to the right and those with small standard deviations to 
the left. This judgment is based on the requirement that at 800 feet 
the curves should have a standard deviation less than 0.3 and at 


deeper depths, they remain below 0.3. With this criteria, artifically 
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Hanie 1, Lewit's (1972) standard deviation (SD) tables for positions 
of interest. Included under total is the number of observa- 
tions that were used to calculate the standard deviation, 
given in degrees centigrade. 

fosition A B C D a 


Total SD S10: Lota 
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Figure 2, Graph of normalized standard deviation for data 
collection points. Values are obtained from Lewit's 


(1972) tables. Figure A is for areas having small 
standard deviation, figure B is for areas with large 
standard deviation. Standard deviation given in a ee 
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imposed by this investigator, points A, B, C and J are judged as 
having a large vertical temperature variability and while at points 
D, E, F, G, H and I, the temperature variability is expected to be 
low. AS expected, points near strong currents, e.g., the Gulf 
Stream, show large variability while points in relatively calm 
waters e.g., the Sargasso Sea, show low variability. 

Table II shows the comparison of the FNWC climatology 
(FNCL) and NAVOCEANO climatology (NAVOCL) with the actual mean 
temperature profile for the three areas A, B, and J. Differences 
between each climatology and the mean for each 12-hour forecasting 
period are given along with the normalized figures. Normalization 
1S accomplished by dividing the difference at any level by the 
maximum difference of either column, 

Normalization is computed to establish a percentage compari- 
son between climatologies of FNWC and NAVOCEANO with respect to the 
mean temperature profile. Figures 3 a,b,c and d are graphs of the 
normalized differences with linear interpolation between points for 
positions B, A, and J. Were a climatology to agree perfectly with 
the mean, the curves would rest exactly on the zero ordinate. The 
aoe Size of the normalized difference at any level is a measure 
of the error associated with that climatology such that 1.0 is the 
Maximum relative error. Thus an area between the curve and the 0.0 
normalized difference line 1S an indication of the error present. 


AREA FNCL - AREA NAVOCL 


= PERCENT DECREASE IN ERROR al 
AREA FNCL = ; ~ a 


Equation (1) is used to give a quantitative measure of the decrease 


in error associated with the FNCL. 
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Table 


Fosition A 


aes 


Comparison of FNCL and NAVOCL with mean temperature pro- 


files. 


Differences are normalized. Number of observations 


upon which the mean is calculated is given in parenthesis. 
Mean profiles are given for indicated Zulu time periods. 
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Figure 3A. A graph of the normalized difference between mean 


temperature profile and NAVOCL and FNCL for position 
B at 0800 local time. 
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Figure 3B. A graph of the normalized difference between mean 
temperature profile and NAVOCL and FNCL for position 
B at 0900 local time, 
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Figure 3C. A graph of the normalized difference between mean 


temperature profile and FNC and NAVOCL for position A. 
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Figure 3D. A graph of the normalized difference between mean 


temperature profile and FNCL and NAVOCL for position J. 
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With this comparison basis, at position A (Figure 3c) FNCL 
yields a 14% smaller error than NAVOCL. Contrarily, for position J 
(Figure 3d), NAVOCL yields a 35% error reduction. For all four 
forecast periods for position B (Figure 3a,b) FNGL shows a substantial 
decrease in error compared to NAVOCL. Specifically, it yields a 
29%, 52%, 51% and 60% decrease in error respectively. 

For position J only two observations are available to obtain 
the mean temperature profile. For the other positions the means are 
based upon a greater number of obServations. Thus for position J, 
that NAVOCL is closer to the mean profile is probably not significant. 
For all remaining forecasts, FNCL had substantially less error than 
NAVOCL. It is concluded from this trend that in areas of expected 
large standard deviation, FNCL is a better approximation to the mean 
temperature profile. 

Analysis of the profiles for positions D, E, G, H and I 
(previously found to be in low standard deviation areas) are listed 
mie bable Ili. As before, the differences between the individual 
climatologies and the mean temperature profile are normalized. These 
normalized values are plotted in Figures 4a-e for each position. 

An area analysis for error is performed as for high standard 
deviation positions. In the low standard deviation areas, FNCL 
showed less error than NAVOCL in every case except position I. 

For position I, NAVOCL had a 38% smaller error than FNCL, the reason 
for which is unknown to this investigator. For positions D, E, G 
and H the results are listed below. Percentages given are the 


percent of decreased error experienced by FNCL over NAVOCL. 
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Table 111. Comparison of FNWC and NAVOCL with mean for the zulu 
time period indicated. Differences are normalized. 
Number of observations upon which the mean is calcu- 
lated is given in parenthesis. 
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Figure 4A. Normalized difference between FNCL and NAVOCL and 


the mean temperature profile for position D. 
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Figure 4B. Normalized differences between FNCL and NAVOCL 


and the mean temperature profile for position E. 
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Figure 4C, Normalized difference between the mean temperature 


profile and the NAVOCL and FNCL for position G. 
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Figure 4D. Normalized difference between mean temperature profile 
and the FNCL ahd NAVOCL for position H. 
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Figure 4E. Normalized difference between mean temperature 
profile‘and the NAVOCL and FNCL for position I, 
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POSITION Bio. SECOND 


FORECAST FORECAST 
PERIOD PERZOR 
D 70% 15% 
E 30% 55% 
G 37% 1% 
H 82% 53% 


As with the high standard deviation areas, FNCL is considered a 
better approximation to the mean temperature profile than NAVOCL 
based upon the strong trend observed. 

2.  ACCustic Predigtions 

Comparisons of temperature profiles and related conclusions 
about their validity is insufficient for a complete analysis. Since 
temperature profiles are not the final product, analyses of the 
acoustic prediction based on each profile are warranted. 

Several methods of comparison of the acoustic prediction 
based on climatology and the mean BT data are available. The para- 
meter used by the fleet for any such analysis is range; therefore, 
this analysis will center on predicted range. The standard set. of 
input data for use with the propagation loss program are: 

wave height = one foot, 
frequency SOTHz; 


source depth = 60 feet, 
receiver depth = 60 feet. 


I 


With this approach all remaining parameters (bottom depth and bottom 
roughness) are held constant except the temperature profile. This 
allows for a comparison of the affects that various separate tem- 
perature profiles have on acoustic propagation. 

The publication, ASW Oceanographic Environmental Services 


(NAVAIR 50-16-24, 1973)explains the method and logic behind the figure 
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of merit (FOM) determination from a propagation loss curve. This 
procedure is followed here wi.th a FOM value of 90 db/ubar assumed 
for all cases. Using the ordinate value of 90 db as a guide, the 
range to the first intersection with the propagation loss curve is 
called the direct path range. Should convergence zone influence 
be strong enough, the 90 db line curve will intersect the propaga- 
tion loss curve at least twice more. The range between the second 
and third intersection is the width of the convergence zone or 
annulus, and the range from the source to the leading edge of the 
annulus is the range to the first convergence zone. Figure 5 shows 
an example of the method for range determination. These three 
values aré@ of tactical importance to the fleet; thus, they are 
pertinent analysis parameters. 

Table IV is constructed for positions B and J, two positions 
in areas of high standard deviation, and H and E, two positions in 
areas of low standard deviation. Ranges for the above three para- 
meters are given. Clark (FNWC, personal communication) indicates 
that due to a computer roundoff of range information, accuracy only 
to the nearest nautical mile can be expected. No greater accuracy 
is dieshaee nor possible by normal computer output means. Such condi- 
tions force the use of a band of ranges vice a single range. Com- 
patibility of either climatology acoustic forecast with the acoustic 
forecast based on the mean temperature profile is achieved should 
any portion of the band of ranges computed from a climatology coin- 
cide with any portion of the band of ranges computed from the mean 


temperature profile. 
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Table JV. Acoustic range analysis for high standard deviation 


positions (B and J) and low standard deviation pro- 
sitions (HW and E). (1) indicates incompatibility and 
(C) indicates compatibility with ranges as computed 
from the mean BT profile. 


Position B FOM = 90 db 
Acoustic Forecast 

Prediction Time 

Based Upon Period 


Direct 
Path (am) 


Acoustic 
Prediction 
Based Upon 


Forecast 
Time 
Period 


Path (nm) 


Forecast 
Time 
Teriod 


Direct 
Path (nm) 





NAVOCL 
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Despite the relative Superiority of FNWC temperature pro- 
files for positions B, H and E, it is necessary to analyze the 
three range predictions separately since different influences 
affected the development of cach. In this way, errors inherent 
in each system can be examined. 

For pesition H, a low standard deviation area, the direct 
path range from NAVOCL is compatible with the mean, but not so for 
the FNCL direct path range. At position E compatibility exists 
at all ranges for both climatologies. The high standard deviaticn 
positions, B and J, show compatible direct path ranges in the 
final forecasts, but are incompatible in the initral torecast, 

i 6. JB EScuss1 on Of Results vA 

One of the main contributers to the magnitude of the 
direct path range is the thickness of the MLD. For cxample, Harvey 
(1972) shows that at 1700 Hz the energy loss for a MLD of 100 feet 
is 2.3 db/nm whereas for a MLD of 200 feet the loss is 1.3 db/nm. 
This MLD thickness effect is significant Since over a typical direct 
path range of 5 nm, 5 db difference occurs. Thus an error in MLD 
thickness of a few tens of feet can appreciably alter a forecast. 
Also, the MLD is a function of space and time. Figure 6 shows the 
MLD from analysis of BT data taken in SEACONEX. The crosshatched 
area depicts the MLD. Both time and range effects are present. 

The cross-section was composed of BT data collected at different 
times and places. The horizontal distance was 32 nm and a time 

of two hours to travel the distance. The variability of the thick- 
ness of the MLD is clearly shown. In places, a change of 50 feet 


in the MLD can occur over a distance of only a single nautical mile. 
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Because of its larger data base, the FNCL was expected to 
show a bettcr correlation to the actual mean MLD than did NAVOCL. 
In general, the temperature profile analysis supports this. Yet, 
Pewee direct pacli Lange 1S a strong function of MLD (Urick, 1967), 
the acoustic forecast indicates no such sSuperiority.; Neither cli- 
matology showed any superiority in acoustic compatibility. For 
positions B, J, and E almost identical records of compatibility 
are achieved. Only for position H did a discrepancy exist. 

Any climatology is, in essence, a mean with an associated 
standard deviation. The actual real-time MLD could lie anywhere 
within certain limits about that mean MLD. The size of the limits 
is dependent upon the magnitude of the standard deviation. For 
example, should the distribution of MLD be Gaussian, there would 
be a 68% probability that the actual MLD would lie within one 
standard deviation of the mean MLD. The actual MLD then will lie 
within an envelope about the mean MLD. 

For a short time period, e.g., one or two days, the actual 
MLD can lie anywhere within that envelope. As long as the FNCL and 
the NAVOCL acoustic predictions lie within that envelope, they have 
equal chance of correlation to the actual MLD. Over a long time 
period the actual MLD will vary about the mean MLD. It is the 
opinion of this investigator that even though a particular climato- 
logy over a lceng time period may provide a better representation of 
the actual MLD than another climatology, over any short time period 
no superiority is observed. Direct ranges predicted by any climato~ 
logy must be considered as a "first guess” until additional data 


LS aVallable. 
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The range to the first convergence zone and width of the 
annulus is a function of the shape of the SVP and the magnitude 
of the gradient of the thermocline. Again it iS expected that the 
larger data base and more sophisticated thermal structure of FNCL 
would allow the ranges associated with it to be closer than the 
ranges associated with NAVOCL when compared with ranges obtained 
by the actual mean temperature profile. As with the direct path 
case, no such superiority was seen in the acoustic range predic- 
tions. For these ranges, the FNCL showed excellent compatibility 
with the mean at all positions except the initial forecast of 
position J. Even for this case, the ranges forecasted were no 
farther from the mean forecast than those of NAVOCL. For the 
NAVOCL, incompatibility 1s experienced at positions H and J. 

Again this investigator considers that conditions on a 
daily basis are too variable for either climatology SVP to yield 
consistently better results than the other. Local effects covld 
influence the mean temperature profile significantly from the 
climatological mean. On the average the results of any correla- 
tion between any two ranges is considered only a temporary effect 
and not to be expected for subsequent reports. This is seen for 
Dost te1on ov; which has no correlation for FNCL in the first fore- 
cast, but exact matching in the second. Figure 6 shows a tongue 
of cold water that extends from depth to near surface. Acccrding 
to Nan-Titi (1967), this tongue may be only a local, short-lived 
effect (a few days) that would be reflected in a local mean but not 
in the Diimatcilectm This is further evidence that a monthly mean 


cannot be expected to correlate exactly with daily observations. 
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This investigator believes that for the purpose of daily 
forecasting no advantage is gained from the sophisticated approach 
to climatology used by FNWC for the case of a transiting ship. 
Should a lengthy time (30 days or more) be spent in a Single area, 
the FNCL forecast may prove to give better results than NAVOCL. 
When only a few days are spent in an area, the daily fluctuations 
about the average will make correlations with any climatology only 


a eirance effect... 


Cee PNOPTIC COMPARISON 
1. FNWC Synoptic Temperature Profile 

Every BT taken during "SEACONEX" was sent to FNWC, where 
it was then included in the mean synoptic temperature profile as 
explained in Chapter I. Propagation loss curves were run for the 
positions of interest (same positions as used in previous climato- 
logy analysis) for each 12-hour analysis period. 

The FNWC thermal analysis forecast at O0O000Z is intended to 
be valid during the period 0000Z-1200Z. This profile was compared 
to the mean of all BT's taken during the interval] 00002-12002. 
This mean, constructed of on-scene BT's, was considered to be the 
actual mean for the period. The same procedure is used for the 
1200Z forecast. No control is placed on the generation of the mean 
profile in that ail BT's are included and.no weighting or biasing 
is used. 

Table V lists the temperature profile as analyzed by FNWC 
for the time periods indicated. The notation in the table 06002 


should be interpreted as 0000Z on the sixth of June. A major 
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Table V. Temperature (C) at various levels for FNWC mean temper- 
ature profile and the actuel mean temperature profile. 
The difference between the two are also given. Zulu 
time periods are given, 


Position A 


Position B 
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Senvesloucion tO this profile,.according to Chapter II, is the BT's 
taken in "SEACONEX." However, BT's taken during the 0000Z-1200Z 
period are not utilized by FNWC until the 1200Z analysis. Thus 

the influence of BT's taken during one forecast period is not felt 
until the next forecast period. Conditions that influenced those 
BT's may not be present during the next forecast period. Hesse 
(1973) suggested that ideally as more BT data are sent to FNWC to 
be included in the mean temperature profile, the analyzed mean 
should better approximate the actual mean temperature profile. On 
the premise that temperature changes in the ocean are slow (Nan=-Titi, 
1967) the BT data sent to FNWC could be used despite the time lapse 
between taking a BT and its utilization. 

One way to test the validity of the assumption that increas- 
ing the BT data base for FNWC analyses will force the temperature 
profile to converge to the actual temperature profile, is to compare 
two actual forecasts for the same position. Table V shows the differ- 
ence between the FNWC mean temperature profile and the actual Mean 
profile. According to the above assumption, this difference should 
decrease with each succeeding forecast. Figure 7 is a plot of the 
difference between the FNWC mean tempcrature profile and the actual 
mean temperature profile. 

Only for position E can a definite pattern be seen for con- 
vergence at all eight levels, i.e., at all levels the difference 
values for the later analysis were smaller than those for the initial 
analysis. For positions B, D, G and H 46% of all observaticns showed 
a greater difference in the later analysis as compared to the preced- 


ing one. Conversely, 45% have a lesser difference and 9% experience 
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Figure 7A. Differences between FNWC mean temperature profile and 
actual mean temperature profile are plotted for the 
time periods indicated at position B, 
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Figure 7B. Difference between FNWC mean temperature profiie and 
actual mean temperature profile are plotted for time 
periods indicated at position D. 
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Figure 7C, Differences between FNWC mean temperature profile and 
actual mean temperature profile are plotted for the 
time periods indicated at position G. 


Difference (°C) 


iO , 220 ; 


¢ 110000 
© 111206 





47 





Figure 7D. Differences between FNWC mean temperature profile and 
actual mean temperature profile are plotted for the 
time periods indicated at position H, 


eee - wa 


100! 


¢ 120000 
O 121200 


400! 


600' 


800! 


1200'& oO : 


48 





Figure 7E. Differences between FNWC mean temperature profile and 
actual mean temperature profile are plotted for the 
time periods indicated at position E, 
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no change. No pattern of genera] divergence or convergence for 
any level can be established, 

Limited data do not allow any definitive conclusions as to 
convergence or divergence. The data that are available show that 
even in areas of small temperature variability, the assumption that 
convergence to the mean temperature profile can be expected with 
increasing data is unfounded. Apparently the problem of time delay 
from observation to data input by FNWC is critical. The temperature 
profile analysis is updated by FNWC every 12 hours. Therefore, 
any BT's taken could be delayed up to 12 hours or more before 
entry into the analysis. Small scale, local anomalies, which may 
have influenced the BT, can change rapidly making the BT informa- 
tion erroneous. This is compounded by the problem of a transiting 
ship Since a spatial change will also come into effect. 

2. FNWC Acoustic Forecast 

Propagation loss profiles for positions H, E and B were used 
to establish a range analysis based on the FNWC synoptic mean tem- 
perature profile and the actual mean temperature profile. The pro- 
eedure was similar to that used in Chapter II. According to Figure 7, 
at ee B the temperature profile converged for the upper 100 feet, 
but was divergent below 100 feet. Position H was an irregular profile 
of convergence and divergence. An entirely convergent pattern was 
experienced at position E. Analysis of these different temperature 
patterns will permit one to study the effects Gna rae torent profiles 
have on the acoustic ranges. 

Table VI shows the forecasted acoustic ranges based upon 


actual mean temperature profiles and FNWC mean profiles. All FNWC 
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Table VI. Range compatibility between acoustic forecasts based 
on a mean synoptic temperature profijJe and the FNWC 
mean temperature profile. (C) denotes compatibility, 
(I) denotes incompatibility. 


Position 5 FOM = 90 db 

Acoustic Forecast Direct ~ First Width 
Prediction Time Path (nm) of 

Based Upon Period Annulus (nm) 


BT Mean 0800 
FNWC Mean 0800 


BT Mean 0812 
FNWC Mean 0812 





Pocition 
Acoustic Forecast Direct Width 
Prediction Time Path (nm) ene 
Based Upon Period Annulus (nm) 
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BT Mean 
FNC Mean 


Position E FOM =_9 
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Prediction Time 
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FNWC Mean 


BT Mean 
FNWC Mean 
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forecasts were based on at least one BT taken in "“SEACONEX." An 
analysis similar to the one for ranyes based on climatclogy in Chapter II 
Section B.2 was used. An envelope of ranges waS considered for 

the three parameters of direct path range to the first convergence 

zone and width of the convergence zone annulus. By definition 
compatibility is achieved should any part of one forecast band 

overlap with the other. 

Position B indicates a considerable error in the direct path 
range for the initial forecast period, but compatibility is achieved 
for the range to the first convergence zone and width of the 
annulus. By the time of the second forecast, compatibility is 
achieved for all ranges. Positions H and E show complete ccmpati- 
bility at all times and ranges. 

These developments setae that regardless of whether tem- 
perature profiles converge or diverge to a mean, once FNWC has 
received several (more than one) BT's from an area the ranges pre- 
dicted are generally compatible with those based on a mean synoptic 
profile. Compatibility cannot be assumed for a single BT input as 
seen for position B which was based upon only a single BT input. 

The feiect path range was incompatible with the mean direct path 

aes All later forecasts were based upon at least two forecasts. in 
Since complete compatibility was observed for these later forecasts, 

it is expected that ranges predicted by ENWC are reliable within 
roundoff error limits, provided at least two BT's hee been input 


into the analysis. 
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ee CAPS [/O ANALYSIS 


Bae TEMPERATURE INPUT 

ICAPS operates on a single BT input upon which the acoustic 
propagation loss profile is generated. In order to establish how 
much reliability can be assigned to a single BT, Table VII is 
presented. For the construction of this table, BT's taken at each 
position are averaged and a variance is calculated for every level. 
Four positions, B, D, E and G, have a variance (07) greater than 
1.0°C for at least one level. Positions H and I have at least one 
level which has a variance greater than 0.5°C. Thus only three of 
nine positions show a variance at all levels less than 0.5°C. 

The above analySis indicates that in all geographic positions 
examined, a high degree of thermal structure variance was preSent. 
Thus at a given location the temperature profile can vary signifi- 
Camely Over a relatively short period of time. A change of 0.5 to 
1.0 degrees centigrade at various levels can alter the MLD and the 
gradient of the thermocline presenting different characteristics for 
saerate proriles: 

{ Figure 6 shows that over a very short distance of the ocean the 
temperature profile can vary considerably. For the profile taken at 
position 1 a 40 foot MLD and a moderately steep thermocline are 
observed. At position 2 the trace shows no MLD and a slightly steeper 
thermocline. The separation between the stations is only two 


nautical miles and 8 minutes in space and time. 
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Table VII. Variances (0%) calculated for all positions based 
on actual BT data taken in “'SEACONEX". 


Position A Position B Position C 
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Effects or internal waves, turbulence, and local anomalies can 
he expected to affect any BT profile. The mean temperature profile 
of an area 1S considered as the best representation of the area, but 
the various influences present can limit any single profile in its 
correlation to that mean. Figures 8a,b and c are constructed to 
show the relationship between a mean profile and the individual BT's 
that were used to establish it. Positions A and H show significant 
differences between BT's with respect to gradient and gradient 
changes. For the most part there is little similarity to the mean 
for any Single profile. For position B, curves (4), (2) and (3) 
do exhibit fairly uniform and reasonable similarity to one another 
and to the mean. Below 400 feet all profiles have a gradient very 
Pieces o70.3 °C/100 feet. Even for this case, curve (1) has a 
MLD of 100 feet while the mean profile does not exhibit any MLD. 

If the mean temperature profile is asSumed to be a reliable esti- 
mate for the true mean profile in an area, this investigator con- 
cludes that a single BT is generally an unreliable estimate of the 


mean temperature profile for an area. 


{ 
JB. ICAPS BCOUSTIC OUTPUT 

As with analysis of climatological and synoptic infomnation, 
temperature analysis alone is inconcluSive. The acoustic informa- 
tionege@eraved by the profiles, being the end product, is the final 
rule against which correlation is measured. 

(ioe VIEL 1S constructed for positions A, B, and H of which 
the first two are positions of high standard deviation and the 


latter of low standard deviation. Using the same analysis for range 





Figure 8A, Temperature profiles for the mean and single BT's 
are given for position H. 
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Temperature profiles for the mean and single BT's 
are given for position B, 
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Figure 8C. Temperature profiles for the mean and single BT's 
are given for position A. 
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Table VIII. Single BT input to ICAPS acoustic ranges versus 
the acoustic ranges for the actual mean temperature 
profile. (C) denotes compatibility, (I) denotes 





incompatibility. 
PocleLom 6 FOM = 90 db 
coustic Forecast “Direct birst Width 
Prediction Time Path (nm) CZ (nm) of 






Based Upon | Period Annulus (nm) 
















0812 = 3) 35 - 37 (C) 4-6 (C) 
0812 1 - 3 35 0= SF 4 - 6 
FOM = 90 db 
Forecast ~ Direct Width 





Time Path (nm) 


of 
Period Annulus (nm) 






Position H FOM = 90 db 

Acoustic Forecast Direct | Perret Width 

Prediction | Time Path (im) CZ (um) of 

Based Upon j; Period Annulus (nm 
1200 i= 3G) 34 - 36 (C) 1 - 3 (C) 
1200 O- 2 34 - 36 eas 





39 





as used previously, it is seen that for the three cases generated, 
complete compatibility for the acoustic ranges 1s experienced. 

" This investigator concludes that at least for the cases shown, 
a single BT can be a reliable indicator for acoustic ranges. The 
variability in temperature profiles has little influence over the 
acoustic transmission loss. / Clark (FNWC, personal communication) 
indicated that low frequency signals tend to be little effected by 
minor changes in the temperature profile. / In other words, minor 
changes in the profile are smoothed out when a low frequency Signal 
is transmitted. Conseguently, the influences of small scale per- 
turbations in the temperature profile will not be realized for low 
frequencies. As long as profiles have the same general shape, 
approximately equal detection ranges can be expected. 

/ With such a small statistical base, conclusions covering any 
large scale area are unfounded. Yet the evidence that is available 
does indicate that for the most part in any single area, a Single 
BT input can furnish the ICAPS operator with reliable range data. 
The ranges calculated can at least be used as a good first approxi- 


Metion te those of the mean within normal error limits. 
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IV. CONCLUSIONS 


A. CLIMA'TOLOGY 

The data presented show that for a large percentage of areas, 
the FNWC climatological temperature profile appeared to be closer 
to the mean temperature profile for the studied areas than the cli- 
matology developed and used by NAVOCEANO for the ICAPS computer 
for the exercise "SEACONEX." It should be remembered that this 
NAVOCEANO climatology was only meant as a first attempt until a 
more sophisticated model could be developed. However, low frequency 
acoustical forecasts based upon either climatology do not indicate 
a clear preference for one system over the other. While local mean 
temperature profiles will smooth, but still incorporate, local 
anomalies, low frequency sound will not respond to small temperature 
PlucCeuatLons. 

This investigator concludes that for a short time period the 
assumptions that FNWC climatology will yield more reliable acoustic 
ranges than NAVOCEANO climatology is unfounded. A transiting naval 
mA can expect nearly the same results from a sophisticated climato- 
logy as well as from a simple one. Forecasts based on climatology 
alone should be used only as an initial approximation and only for 
information purposes. The probability that the predicted ranges will 


Petlecesoctvael Conacitions are uncertain. 


Bs or NOPTIC 
j 
Y FNWC synoptic temperature profiles show no overall evidence of 


convergence to a mean temperature profile. Yet, within limits of 


. 61 





roundoff error, the ranges forecast by FNWC based upon at least 
two BT's proved to be similar to that based upon the mean BT pro- 
file. It is concluded that FNWC synoptic forecasts do indeed 
PLroviae a reliable tactical forecast . 

The single BT input for ICAPS has proven to be a reliable 
forecast agent. However, the user must attempt to calculate a 
mean temperature profile for his area as soon as possible. Until 
a mean profile is established, a single BT input may be used with. 


moderate reliability. 


C. OPERATIONAL CONSIDERATIONS 

Based on information from "SEACONEX," the large time delay (up 
to 12 hours) between the time a'BT was taken and incorporation in 
the temperature profile at FNWC and the time delay (average of two 
hours) between the time the FNWC forecast was made and receipt of 
the forecast, and the fact that in any combat Situation, communica- 
tions with FNWC will be minimal, the dependence upon FNWC to provide 
real-time forecasts to the fleet seems unrealistic. ICAPS utilizes 
the same acoustic models and inputs without the communications 
problems with a very short time delay. It is concluded that the 
ICAPS system is the appropriate real-time system, but FNWC can he 
used aS an accurate backup system. 

For tactical planning purposes, FNWC climatology should be used 
as the source of acoustical range information. Their large choice 
of data presentation techniques and ability to widely vary input para- 
meters make them well suited to laying out the expected acoustic 


Banges rOm any area 1n a style best suited to the situation. 
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Greater detail and lesser roundoff error can be obtained by a 
proper choice of display. At present this option is not available 
to ICAPS. Also, this study showed that for short time periods, 

the superiority of any climatology versus another is unfounded. 

Yet for any longer period of consideration, the much more sophisti- 
cated climatology of FNWC as opposed to the interim climatology of 


NAVOCEANO should present more representative information. 
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Ve RECOMMENDATIONS 


A. STANDARD DEVIATION FIELDS 

In order to supply the ICAPS operator with a scheme upon which 
to base reliability of his forecasts, it 1s suggested that NAVOCEANO 
compile Lewit's (1971) tables in such a way that the expected 
standard deviation from the mean temperature field for each acousti- 
cal area can be made operationally useful. A color code for large, 
moderate or small standard deviation fields is a simple but useful 
guide. Such a presentation will afford the operator an insight into 
how much variability can be expected in his area. Also, he could 
plan his data collection based on such a system, e.g. more BT's 
could be taken in areas of large standard deviation and fewer in 


areas of small standard deviation. 


Be. TRAINING OF ICAPS OPERATORS 

NAVOCEANO and the ASW schools should conduct a short course on 
the proper operation of the ICAPS system. Ships equipped with ICAPS 
should be assigned trained oceanographers with a sense of apprecia- 
tion for real-time problems and the art of good data collection and 


forecasting methods. 


Cs ALTERATIONS TO ICAPS SYSTEM 
BT input should be recorded on a cassette tape for easy handling 
and recovery. Presently a profile must be hand digitized and entered 


into the computer. Also the program should be altered to allow for 
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the input of more than a Single set of parameters at a time. The 
program could solve and display each set separately, which would 


Save consiaqerable time. 


D. FNWC PROGRAM ALTERATIONS 

Presently FNWC calculates a variable constant which is an indi- 
cation of the rapidity with which the mean tempcrature profile may 
vary from its initial state (Clark, personal communication). This 


constant iS based upon the data density (number of observations 


used to form the statistics), the sonic layer depth and the magnitude 


of the sound velocity gradient. The number is then simply sent to 
the user. 

Two basic problems are noticed in this approach, Parameters 
such as expected standard deviation of the temperature profile, and 
actual roughness and slope of the bottom spatial and temporal varia- 
tions in received BT profiles are not includcd in the program. More 
importantly, the user does not know which parameter exerted the major 
influence in the Se eek of the constant. For the first case, 
Ee program should be rewritten to inciude the three parameters 
listed since they exert a large effect on the variability of the 
acoustic profile. Secondly, for the latter problem, the user should 
be informed as to which parameter had the greatest influence on the 
forecast in order to enable him to better deal with the problem. A 


Simple code could be established showing the degree of relative in- 


fluence each had on the acoustic profile. In this way the usefulness 


of a single BT is enhanced in that a user would be able to judge as 
to the variability of the water column, spatial variation, and other 


parameters separately. 
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Therefore, it is the opinion of this investigator that the 
program be rewritten to include additional parameters such as 
those discussed above. Additionally, some method to inform the 


user of the most influential parameter should be developed. 
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APPENDIX A 


FNWC SEA SURFACE TEMPERATURE ANALYSIS 


The basis for the Fleet Numerical Weather Central's (FNWC) sea 
surface temperature analysis scheme is the Fields by Information 
Blending (FIB) program. Figure 9 shows a flow diagram for the FIB 
program. The analysis scheme requires a first guess temperature 
field, ver and a base field temperature, T Figure 10 shows the 
area module for the program. Gradients are given aS wu and Vv 
with associated weights B* and C* (Holl et al, 1971). 

The basis for T. 1S the preceding analysis, T , where 4 
is the analysis period, adjusted slightly toward climatology. The 
climatology, Tos is computed as a weighted mean for each analysis. 
The weight assigned to the current month is one~half,and the preced- 
ing and following months have a combined weight of one-half, pro- 
raced by the hour of the month. Thus the equation for the first 


guess temperature is as follows: 


Tl. = (1-F) *Tr rae! shoal 


G c 


where F is proportional weight assigned to climatology (Holl et al, 
OTL). 
In order to nullify the small scale anomalies, often caused by 
sparse data, a climatic smoothing is accomplished: 
= — le + rt 
B G To) sMoorHED To 


The first guess difference ficla, J Bee 1S given by 


7 ool Toe (HOll-ot al, 1971) 
O B 
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Figure 9, Flow diagram for SST analysis 
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Feeure 10, Grid diagram for SST analysis 
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Observations (subscript n) are transformed into difference 
values: 
Dn ae ene ol ect al, 1971} 

Data preparation 1S a Screening process designed to: (1) elimi- 
nate all reports not containing an SST value; (2) eliminate duplicate 
reports; (3) reject reports outside analysis region; (4) assign 
report reliability, Avs (a function of the contributing variances) ; 
(5) perform an initial gross error theck by analysis the magnitude. 
of os ; and (6) reduce reports weights for quasi-duplication, that 
is, for non-identical reports from the same reporting platform 
(Holt et al, 1971). 

During the assembly phase, each report is assembled at the nearest 
grid point. The report difference value, Cae and its weight, Aas 


are added to the accrued information at the grid point: 


Ne Ni eal 


ae oe Bee oN-a oN 


N refers to the stage of the assembly and no to the single report 
Gielimee al, 1971). 


The blending equation reduces to 


J 


ih % 
= + + 4 
Wen Som *2,m,n72,m,n * 'Bem4, mel Co harm Bo m-1 &mm-1 
#Cy a m%t-1,m! 
where S is defined by > n P x 
= + + + 
Bn 9 sn en ai arn Po m1 en 


A* is solved for using the following equation: 


a 
5m 


* 


_ =] 
Ag pm 7 Ae jm Ag m/ A 
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Heeentesteetie L£1¢St pass through the Assembly, J* Blending, 
and A* Solution phases, a gross error check is now made on each 
SST observation. All reports are re-evaluated considering the new 
J* field. Once again the accepted reports are assembled at the 
nearest grid point and the J* Blending and A* Solution phases are 
completed. This time the output is combined with the Base Field 
Temperature (T,.) to provide the analyzed sea surface temperature 


Prelas (Hol? et al, 1971). 


on 





The 


APPENDIX B 


FNWC MIXED LAYER DEPTH ANALYSIS 


FNWC mixed layer depth (MLD) analysis scheme is descrihed 


by Hesse (1973). The MLD is defined as the lower boundary of the 


turbulent, mixed surface layer or the upper boundary of the thermo-~ 


cline, 


The potential mixed layer depth (POTMLD) is defined as the 


Gepth to which the turbulent mixed layer depth extends. It repre- 


sents the mean MLD about which short-period fluctuations can be 


Expeeced. 
The model has three steps: 
1. Adjusting the first guess--the 12-hour-old MLD analysis-- 
climatology to form Guess l. 
2. Computing the effects of wave action on the MLD, and 
adjusting Guess 1 accordingly. 
3. Adjusting the guess field (Guess 1) to available BT reports 


taken in last 72 hours, 


In the actual program computations the analysis proceeds as 


follows: 


ibe 


Guess 12-hour-old POTMLD 


0 


Guess, = Guess, fot 7 ro. (CLIMAT = Guess), where 


if day = 14, CLIMAT = (0.5 - day/30) + (day/30 + 0.5), 


1f day = 15, CLIMAT = CLIMAT for current month, 


Le L6 day = 30, CLIMAT = (1.5 ~ day/30) (CLIMAT current 
month + (day/30 - 0.5) (CLIMAT next month), 


if day = 31, CLIMAT = same as for 30. 
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Compute MLD; 


MLD = GH =“Osissir- fT 


LE Sol 
waves 600)?! * 


CH * J1 if SST = T 600 


where CH = combined wave height, 


SST = sea Surface temperature, 


= 
Il 


temperature at 600 feet. 


Smooth the guess field with a special smoother. 


Force the guess field to be at least as deep as the MLD 


waves field: 


Guess. = Guess if Guess. 2 MED 
2 ee i mes waves 
= D ps LEwGDeSSs=4, MED, 
waves I: 


Adjust Guess, LLié€ld oO Bie re vee. 


Limit final analysis field by forcing POTMLD 
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Figure 1]A. Basic flow chart for sub-surface temperature analysis. 
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Figure 11B, Basis flow chart for sub-surface temperature analysis 
(continued ) 
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APPENDIX C 


FNWC SUB-SURFACE TEMPERATURE ANALYSIS 


The FNWC sub-surface temperature analysis is run at the 100, 
200-7300; 400, 600, S00%and 1200 feet levels. The 12=hour old 
analysis iS used as a starting point for each level. This is 
adjusted 1/16 of the way toward climatology. After this adjust- 
ment, the figures are further adjusted to the MLD and SST analy- 
sis to produce the guess field. The MLD and SST analysis serve 
as strong anchor points for the subsequent sub-surface analysis. 
(Hesse, 1973) 

BT reports for the preceding 72-hour period are included in 
the analysis, and they are weighted by age on al to 4 scale. 

At each level the report is subjected to a gross error check (GEC) 
by requiring that the report fall within limiting values of the 
guess field. A similar GEC is run on reported gradients by com- 
parison with the guess field gradients. (Hesse, 1973) 

The guess etiatis) we adjusted to all reports which passed the 
GEC. The vertical gradient between the first analysis and the 
final eee from the level above is then slightly smoothed. 
Next, consistency is checked by allowing a positive temperature 
gradient to exist only at selected grid points. Finally, the tem- 
perature range at all levels is restricted to the range of -2 
degrees to 35 degrees centigrade. Figure 10 shows the basic flow 


chart for the sub-surface analysis. (Hesse, 1973) 
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